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We study the response of the nucleon, as a system of three bound (constituent) non relativistic 
quarks, to external (quasi static) electric and magnetic fields. The approach, based on a sum 
rule technique, is applied to a large class of two and three-body interquark potentials. Lower and 
upper bounds to the electric polarizability and para-magnetic susceptibility are explicitly calculated 
within a large variety of constituent models and their values related to the features of the interquark 
qq", interaction picture. The role of three-body forces is investigated in details as well as the effects 

of SU(6) breaking terms in the potential model. Our results can be used to extract the mesonic 
' contributions to the static polarizabilty and susceptibility. The quark degrees of freedom give a 

I , quite sizeable contributions to both and the meson cloud accounts roughly for 30% and 60% of the 

electric proton and neutron polarizability respectively. The quark contribution to the paramagnetic 
^"t ' susceptibility is even higher and the mesonic effects are rather uncertain. 

>—i '. 

<N ■ 

C\| , I. INTRODUCTION 

The electric and magnetic polarizabilities, labeled ole and (3m respectively, are fundamental observables to under- 

■ stand the intrisic structure of the nucleons. They caracterize the ability of the constituents to rearrange in response 
I/* - ) , to (quasi) static external electric and magnetic fields and are as fundamental as other parameters like charge radii 

■ and magnetic moments. 

For the proton they can be measured by means of Compton scattering experiments since the Thomson scattering 
amplitude X/j = — e-e' e 2 /Mjy is modified, for intermediate values of the photon energies (50 MeV < to < 100 MeV), by 
first order corrections due to the nucleon polarizabilities yielding to the low-energy Compton scattering cross section 



-3': ^> = ^>-4(^> 

s : 
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^ E + — (1 + COSf?) 2 + ^LJk (1 _ cos of 



(1) 



\ where da B j d£l is the Born cross section for a proton with an anomalous magnetic moment and no other structure (e.g. 
ref. Q), and e 2 is the fine structure constant, cxe and (3m are the so called dynamic (or Compton) polarizabilities 
containing (within a non- relativistic approach) retardation effects as well as diamagnetic contributions We will 
discuss such effects later on. 

Eq. (Q) shows that forward and backward cross sections are dominated by a E + (3 E and a E — (3 M respectively, while 
experiments at 90° are sensitive to a E only and the polarizabilities are obtained by measuring the deviations of the 
cross section from the Born values. In addition a dispersion sum rule [|J constrains the sum of a E and [3 M : 



W M = ^\ ^#^« 14.2, 0.5 10 ! I'm : . (2; 



where a\ ot (ixi) is the total proton photoabsorption cross section and the numerical value is obtained extrapolating tha 
available experimental data to infinite energy 

A global analysis of the existing experimental data has been recently performed yielding to 

a\ = (12.1 ±0.8 ±0.5) 10" 4 fm 3 ; 0^ = (2.1 ± 0.8 ± 0.5) 10~ 4 fm 3 ; (3) 
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The neutron has no charge and Compton interference experiments to measure a E are not possible^]. Recent 
measurements (at ORNL) have been realized by means of a precise n-Pb scattering experiment 0. The electric 
dipole moment, induced by the nuclear charge on the moving neutron, produces a secon order effect resulting in 
a 1/r interaction proportional to the neutron electric polarizability which can be measured looking at the energy 
dependence of the neutron- 208 Pb scattering cross section. The result is 

a% = (12.3 ± 1.5 ±2.0) l(T 4 fm 3 ; = (3.5 ± 1.6 ± 2.0) KT 4 fm 3 ; (4) 

where [3^ is extracted by assuming a E + = (15.8 ± 0.5) 10~ 4 fm 3 according to the sum rule (§), and the 
retardation effects are included following a suggestion due to L'vov and Petrun'kin. The direct experimental result 
of the n-Pb scattering is a E rather than since the measurement has a static character, and the value is a E = 
(12.0 ± 1.5 ± 2.0) 10 _4 fm 3 which is quite close the Compton polarizability because the corrections for the neutron are 
small. 

For our study we will use non-relativistic approches and, in this case, the corrections assume a rather simple form. 
In fact the Compton polarizabilities can be written (e.g. Friar in ref. |l|]) 

^ = ^<°iL e ^ i°) + 2 L En _ Eo = 

i=l n>0 " U 

= Aa E + a E (5) 



P. 



\(n\^\0}\' 



\i=l / i=l n>0 n ' 

A/3M+/3M| dia + /3Af| para • (6) 



and 

^ (0| (t « J\ 10) - |<0| t ef rf M|0) + 2 £ ^ 

\i=l / i=l n>0 n 

= A, 
In Eq. (|) 

i=l 

is the nucleon electric dipole operator in the z-dircction, i.e. in the direction of the external electric field. The 
coordinate of the charges z[ refere to the center of mass of the nucleon (the relevance of the motion of the center of 
mass in dipole excitations is one of the motivation for using non-relativistic model where the spurious contributions 
can be separated exactly) and AaE is a correction due to retardation effects Aa^ = (r 2 ) c h/3MN and is related to 
the charge mean square radius of the system. 
In Eq. |) 

1=1 1=1 l 

is the nucleon dipole magnetic moment operator, Oi the Pauli matrices. (3m contains both the diamagnetic and the 
paramagnetic susceptibilities and, in addition, a retardation correction that can be written 

A/3 M = -^<0|D^|0>; (9) 

Eqs.(^[^) are valid both for neutrons and protons. 

From a theoretical point of view the nucleon polarizabilities received much attention and, in the last few years, 
many approaches have been developed within complementary frameworks like constituent quark models ]9|-[Tl|] , MIT 
bag models (l2| and its chiral extensions , soliton models Jl4],[l5| , chiral perturbation theory |l6|,[n]] and dispersion 
relation method fH. 



1 Compton scattering on the neutron can be realized by means of deuteron targets, e.g. ref. 
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In the present paper we want to reconsider the non-relativistic constituent quark model which is often believed to 
be inadequate to describe the static response of the nucleon to external fields (cfr. e.g. the discussion in refs. 
As a matter of fact the constituent quark model does not incorporate meson degrees of freedom, and consequentely 
the effects of the meson cloud cannot be included. Therefore one cannot expect to reproduce the exprimcntal value 
of the electric polarizability considering quark degrees of freedom only, since the meson cloud should be responsablc 
for a sizeable part of the electromagnetic reponse. However the exact amount of the quark contribution has not 
been calculated by using potential models which are fitted on the baryonic spectrum, but simply estimated within 
naive models. In the following we develop a sum rule approach to the nucleon polarizabilities which is quite general 
to include two-body and three-body forces between the constituent quarks. In particular three-body forces have 
been recently considered in the context of quark models to refine the baryon mass spectrum predictions and, in 
particular, to reproduce the position of the Roper resonance. Moreover three-body forces are strictly related to the 
gluon-gluon interaction which is one of the fundamental aspects of QCD. In addition we consider consistently the 
retardation effects arising from a non-relativistic approach to the electric and magnetic response. Our results should 
be, therefore, considered as a precise estimate of the quark contribution to the electric and magnetic response. 

The paper is organized in the following way: in section [n| we describe the sum rule approach to the linear response 
of a non relativistic system and introduce a technique to evaluate various upper and lower bounds to the polarizability 
sum. The relation between the interquark potential model and the electric and magnetic sum rules is discussed in 



detail in section III where the Isgur-Karl model, three-body hyperradial and two-body plus three-body potential 
models are considered and discussed for both electric and magnetic excitations. The Hamiltonian is solved exactly by 
means of the Schrodinger equation and the wave functions of the nucleon used to calculate the sum rules. The role 



of the confining, hyperfme and three-body forces is investigated in detail and the results summarized in section IV 
where the comparison with experimental data is also discussed. 

II. SUM RULE APPROACH TO NUCLEON POLARIZABILITY 

A quantity of the type 



m -i( e ) - L -e—^ ( 10 ) 

n>0 " U 

is involved in both the expressions (j^) and (^|) for the electric polarizability and the (para) magnetic susceptibility 
(0 corresponds to the electric dipole (fj]) and the magnetic dipole operator (||) respectively). A direct evaluation of 
to_i(0) from Eq. ( JTo| ) involves all the complications of the excitation spectrum of the system (energies and wave 
functions) and it cannot be easily estimated truncating the excitation space fl9fl . However one can construct upper 
and lower bounds to m_i [0,pllJ2(T] by using the positivity of the corresponding strength distribution 
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5eH = ^|(7i|e|0)| 2 ( S(^-(^ l -So)) , (11) 
and few (positive) moments of such distribution (p > 0, integer) 



m p (Q) = / du uj p Sq(ui) = 
Jo 

= £|(n|9|0)| 2 {E n -E a f = 

n>0 

= (o\eHH -E re\o)~5 p0 \(o\e\o)\ 2 . (12) 

The last expression has been obtained by using the closure property of the eigenstate \n) and can be expressed in a 
simple form for the first few moments, or sum rules ( for reviews on the sum rule techniques see e.g. refs. p2|| ) 

mo(9) = (0|ete|0)-|(0|9|0)| 2 , (13) 

mi (Q) = |<o| [et,[H 0) e]]|o) , (14) 

ma(0) = i<0|{[et,iTo],[flb,e]}|0) , (15) 

"*s(e) = ^<0| [[et,iT ] ,[H ,[H ,Q]]] |0) , (16) 
etc. . 
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The sum rules (13 - 16) involve the ground state wave function only, and commutators (anticommutators) of the 
nucleon Hamiltonian (Hq) with the excitation operator (0 = D' z , fi z ). In addition the dynamical aspects embodied in 
the Hamiltonian enter in a more and more complex way increasing the order of the sum rules and the explicit calculation 
are of increasing complexity. In the next sections we show how one can construct bounds to the polarizability by 
means of sum rules of different order. 

A. Lower bounds 

In the limiting case of a strength distribution concentrated in a very narrow region (<5-like distribution), the polar- 
izability can be easily evaluated considering sum rules of positive order only. In fact 

,2 



m -i( ) = ~^77^ = \/mi(e)m- 3 (e) etc. (17) 

However the strength is generally distributed over many eigenstates and its spread affects the lower and upper moments 
in a different way, and the equalities ( |l7|) become rather inequalities 



m_x(9) > ; m_i(0) < ^(9) m_ 3 (6) etc . (18) 

An elegant way of taking partially into account the effects of the spreading and the width of the strength distribution 
has been proposed by Dalfovo and Stringari [^(j) minimizing, with respect to the parameters a and b, the inequality 

doJ^^- (l + au + buj 2 ) 2 >0 (p integer) , (19) 

based on the positivity of Sq(lu) [j. One obtains 

m 2 „ , , (0) 1 

m_ p (9)> - p+1 n i A(n , /r(n , (20) 
m_ p+2 (0) 1 - A p (9)/r p (9) 

where 

A p( 0)=f^-^y , (2D 



and 



r P (9) 



m- p+ 2 V m -P+1/ w-p+l 
A lower bound to the polarizabilty sum rule is obtained for p — 1 



(22) 



m- X 0>-^-( - — — — — . 23 
m 1 (9) 1 - Ai(9)/ri(9) 

In addition to the previous relations, other bounds can be obtained by using Schwartz inequalities: an example 
involving the first odd sum rules is 

m-m > ^§ (24) 
m 3 (0) 



based on the conditions 



m p (0) < m p+2 (0) 



m p _ 2 (0) m p (0) 

valid for the integer values of p = , ±1 , ±2 .... etc. Eq. (pj) is obtained for p = 1 



2 The authors of ref. pOl considered the case p = 1 only. 
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B. Upper bounds 



Upper bounds to the polarizability are easily found considering that the energy transfer lu cannot be smaller than 
the energy difference between the first excited state of the system and the ground state, i.e. li > E\ — Eq (= Eiq). 
The simplest way is just using the definition (vim and applying closure 



m rei - V K"IQ|Q>I 2 < y l(»ie|Q)l 2 _ ™o(Q) 

U 1 4- ' E n -E ~ ^ n E 10 E w 

n>0 n>0 



Corrections to the previous bound can be considered following again ref. |20| 



f 

Jo 



aid (1+7CJJ < 



which leads to 



m_i(6) < 



m (9) 
Eio 

m (9) 
-Bin 



(1 + 7CJ) 



in 



(26) 



(27) 



1-^ ^-£ 10 I^-Exo 

ml \mo / \ TO i 



A(9) 



(28) 



An even more stringent upper bound can be obtained again from Eq.([19|): in fact, for p = 2 one obtains 

m 1 (e)m (6) 



m_i(6) < 



m 2 (6) 
mi mo 



1 



m n m 2 
mi 



- 1 



/ m 2 m_ 2 
V ml 



- 1 



m 2 



S(6) 



(29) 



Note that the root argument is always positive (due to Scwartz inequalities and the relation (|l^) for p = 1). Eq. ( |29| ) 
involves inverse quadratic energy weighted sum rules. 

Finally, from the Schwartz inequalities (Eq), one can obtain also an upper limit for m_i from 



rn-i(e) < Vmi(e) m_ 3 (9) 



(30) 



which involves the inverse cubic moment. 

m_2(9) and m_3(0) have no closed form in terms of commutators and/or anticommutators; due to the quadratic 
and cubic energy power at the denominator they can be easily estimated including the first excited states only (in 
the following we will make use of the first two states). 



III. SUM RULES AND INTERQUARK POTENTIAL MODELS 

The main ingredients entering the sum rules are: i) the nucleon ground state, and ii) the corresponding Hamiltonian 



H =T + vW q + v}%_ q , 



(31) 



which includes, in principle, two and three-body interaction terms. 



A. Harmonic oscillator potential model 



The most simple choice is the harmonic oscillator two-body confining potential, ft has analytic solutions and 
supplies a convenient classification scheme of the baryon resonances. If we assume a (spin-independent) harmonic 

form (V q ^} q _ q = 0), V q _ q — 1/2 K J2i<j( r i~ r j) 2 = 1/6 mu)\ J2i<j( r i~ r j) 2 > ^ i s rather easy to derive the sum rules 
for the electric transitions by noting that all the excited electric dipole states are degenerates at E D = E® = 1 ftwh.o. 
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and that 
form El 



D' V {2) 

^ zi v q—q 



V, 



(3) 



O.As a consequence the linear energy-weighted sum rule assume a model independent 



D't,[H ,D' z ] |0> = ±<0| D'l,[T,D' z ] |0> = 



and the sum rules can be written 

m p (D' z )=m 1 (D' z ) (E^ 1 



3 m 



;2 

3m 



-<o. 



e a 



2p-2 



m )J 



P = 0,±l,±2,etc. 



(32) 



(33) 



with aj; = mwh.o.- When p = — 1 one derives the prediction of the harmonic oscillator model for the electric 
polarizability of neutron and proton 110] 



n.P _ 



2m- 1 (D' z ) = -e' 



M 



N 



'h.o. 



(34) 



The previous equation can be expressed in a more general form as function of the proton charge mean square radius 
(r£)ch = l/aL.> 



■e 2 M N (r 2 p ) 2 ch , 



(35) 



which can be derived within a variational approach to the problem as it will be extensively discussed in section IV A . 

The paramagnetic susceptibility vanishes in the h.o. model since the energy variation to flip a quark spin is zero 
(the interquark potential dose not depends on the spin degrees of freedom and therefore the A's mass does not differ 
from the nucleon mass). From a sum rule point of view we note that, since [Hq, /i z ] = in the h.o. potential model, 

TO i(/ i z) = TO 2(Mz) = m 3(Mz) = Oj while mo(/i z ) = | (5^) (cfr. section [II C 2) and the bounds previously introduced 
are not well defined. 



B. Potential models and variational approaches 



In more sophisticate versions of the interquark interaction, the confining potential consists of the h.o. part plus 
terms which removes its degeneracy. This is the case of the Isgur-Karl (IK) model which contains an unknown 
[/-term added to shift the energies of some states. The model includes a delta-like hyperfine interaction derived 
from the one-gluon-exchange gg-potential. The baryonic states (obtained diagonalizing the Hamiltonian within a h.o. 
basis) result in a superposition of different SU (6) configurations and lead to a rather good description of the baryonic 
spectrum. 

A sum rule approach to the nucleon polarizabilities can be developed also for that kind of potential and we discussed 
results on in a recent work |ll[] . However the use of sum rule technique is limited to potentials and wave functions 
which are well defined and self-consistent. The IK model, making use of an unknown potential term, prevents explicit 
calculations of the sum rules. As a consequence the approach proposed in ref. Jll| was limited to few sums (namely 
mi, 7713 and m-3) and to the bounds 

2 < a E < 2 y/ mi (D' z )m- 3 (D' z ) , (36) 

because they are the only ones which still contain some informations. In fact mi is model independent and the 
unknown (central) [/-potential is irrilevant. m_3 can be evaluated explicitly making use of the first two excited 
electric dipole states of the nucleon spectrum (Sn and -D13) and therefore the upper limit is well defined. The lower 
bound, even if it is written as a combination of sum rules, defines a variational approach to which assumes that the 
global effect of the external electric field (£) is a rigid translation of the u-quark wave function against the c£-quarks 
without any additional deformation |p5| . As a matter of fact variational and sum rule approaches are often related 
and we discuss this point in some detail. 
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1. sum rules and wave function deformations 



Let us assume that the wave function of the polarized nucleon can be obtained from the unperturbed ground state 
by means of a unitary transformation shifting charges in opposite directions in a rigid way, one can write a set of trial 
states of the form 

\<t>e) = e^ A \0) , (37) 
A is an anti-hcrmitian operator and r\{£ ) the variational parameter. Assuming A — [T, D' z ] one has 

\<t> £ ) = e " (£) [ T '^] |0) = e -^( f ) 2m E i e 'fM | ) . (38) 

The states \4> £ ), driven by the translation operator p' z are shifted, in the z-direction, in a rigid way. For weak external 
fields rj(£) — > and one can expand \<j) £ ) to evaluate the variation of the total energy. One gets 

Etot = (cl> £ \H - D' z £\<j> £ ) 

- (0|tfo|Q) - y (0| [A, [H 0l A]} |0) - £ V (0\ [D' Z ,A] |0) (39) 
and the expectation value of the induced dipole defines the polarizability 

(<t> £ \D' z \ci> £ )-(0\D' z \0)=r ] (a\[D' z ,A}\a) = a E £ . (40) 
Minimizing the energy variation ( p9| ) one can extract the equilibrium value of the variational parameter 

1 (0| [A, [H , A]] |0) 1 ) 

and substituting this expression in Eq. ( flof ) one finally obtains 

((0|[^]|0)) 2 



OLE 



<0|L4, [H o ,A]]\0) 

((0\[D'AT,D' z }]\0)) 2 
)| [[D' z , T], [H , [T, D' z }}} |0) 



(42) 



Eq. ((42|) represents a lower bound to the polarizability |26) since the solutions of the variational equation for \cf) £ ) 
are restricted to the wave functions obtained just shifting the ground state charge densities of the u and d quarks in 
opposite directions as described by Eqs. (^) and (^). However when the Hamiltonian Hg of the system in study 
has the property [Ho,D' z ] — [T, D' z ] because the potential commutes with the electric dipole operator, Eq. ( (42| ) can 
be written as combination of sum rules 

((0|[^,[ff ,^]]|0)) 2 _ ? rnl(D' z ) 
E mD' z ,H ],[H ,[H ,D> z }m " rn 3 (d' z ) ! 1 } 

which is valid for all the models without velocity dependent (or non-local) potentials, and is equivalent to the lower 
bound disussed in section II A and used in ref. . 

Other possible variational approaches have a close relation with sum rules and we discuss a second example which 
assumes that the wave function of the system, under the influence of the electric field, can be approximated by means 
of a simple deformation driven by an operator F such that 

|#) = ^= [1 + aF] |0) . (44) 

a is a variational parameter, |0) the unperturbed ground state of the nucleon, and F = J2i f( r 'i) a single particle 
(local) operator depending on the (intrinsic) quark coordinates only. y/N £ ia a normalization factor which ensures the 
condition (ip£\ip £ ) = 1 and /(r'i) a function of r'i to be guessed in order to introduce physical "dipole" deformations 
on the wave function. We note that F commutes with the dipole operator so that [F, D' z ] — 0, (0|F|0) — because of 
parity and the "dipole" character of the function F, and N £ = 1 + a 2 (0|F 2 |0). 
Calculating the total energy variation one gets 
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-» (l -a 2 (0|F 2 |0)) «0|JEf |0) + a 2 (0|Fff F|0) - 2af (0|F£>;|0)) 



= #o + y<0|[F,[iT ,F]]|0>-2a£<0|FD;|0) , (45) 



and 



<<M^I<fe) - (OPHO) = 2a(0|F£>;|0) = a E £ . (46) 
Minimizing ( |45| ) and substituting in ( (fuj ) the obtained equilibrium value of a, one has 

ag = 4 mFD'Mf (47) 



The previous expression is similar to Eqs. (|42j), ([43|) and can be evaluated by using canonical commutation relations^. 
A quite natural and simple choice for F is F — D' z . In this case 

((OI^^IO)) 2 _ ml(D' z ) 
E (0|[U;,[H ,^]]|0) mi(D' z ) ' 1 > 



and a simplified lower bound is found. It is rather similar to the bound (|23|) and can be obtained from the constrain 
© in the limit p = 1 and b = 0. 



C. Hyperradial potentials and three-body forces 

Just at the opposite side of the simple two-body harmonic potential, one can locate the three-body force model 
(TBM) recently proposed by Ferraris, Giannini, Pizzo, Santopinto and Tiator ^7|. The model assumes that the 
interquark potential can be written as an hypercentral potential, i.e. a potential depending on the hyperradius 
only. The idea of hyperradius (£) is introduced in the so called hyperspherical formalism |28| together with the 
hyperangle (<fi), to define the hyperspherical set of coordinates in a six-dimensional space: Cl p , Q\, £ and <j), with 
£ = \J p 1 + A 2 , = arctan (£) and p and A are the absolute values of the Jacobi coordinates p = (v\ — Y2)/\/2, and 

A=( ri +r 2 -2r 3 )/V6. 

Under the assumption of an hypercentral potential the intrinsic Hamiltonian Hq can be written 

H ° - -L w + v » + 1 " 3, «» —t + U + ^) + v,3 '«» ■ {m> 

where L 2 {Vl)/^ 2 is the analogous, in six dimensions, of the three-dimensional centrifugal barrier and Q, embodies 
Vlp,Q,\,(j). The eigenfunctions of the grand-angular operator L 2 are called hyperspherical harmonics and denoted 
by YLi(Q). They are written as products of spherical harmonics in Q p and £l\ with angular momentum l p and l\ 

(corresponding to the coordinates p and A) and Jacobi polynomials in the hyperangle cf> p8| . They form a complete 
orthogonal basis in the space of the five-dimensional functions of Sl p , Sl\ and cf>. The eigenvalues of L 2 are —7 (7 + 4), 7 
are called the grand-angular quantum numbers and are given by 7 = 2 k + l p + 1\ with k integer and non-negative. As 
long as the potential is hypercentral the complete wave function ^(^Jl) can separated as ^(^, fi) = il} Vj y(^)Y\~\(d) 
and the hyperradial wave function Vv,7(£) i s solution of a Scrodinger equation. For fixed value of 7, the label v 
indicates the number (y + 1) of nodes in the wave function. 

The hypercentral character of the potential means (in general) that the interquark interaction has a genuine three- 
body character, in the sense that the coordinates of a specific pair cannot be disentangled from the third one [p9| . 
The idea of a three-quark force is related to the non-abelian character of QCD, in particular to the existence of a 
direct gluon-gluon interaction which represents a fundamental justification for the introduction of three-body forces 



3 The two example given in this section could be unified within the unitary transformation formalism of Eq. ([^), consid- 
ering the operator A as a many-body operator satisfying the relations (0| [A, D' z ] |0) = (0|-FD^|0} and (0| [A, [Ho, A]] |0) = 
(0\[F,[H ,F]}\0). 
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. From a more phenomenological point of view it has been argued [ pl|]32] ] that the contribution of a three-body 
force is of some help in solving the "Roper puzzle" |33|] giving a consistent explanation to its relative position in the 
spectrum. 

For a given hypercentral potential the eigenvalue equation has to be solved numerically obtaining energies and 
hypercentral wave functions which, depending on the hyperradius only, result to be completely symmetric. The 
correct global symmetry of the state is obtained combining ip Vj -y(£) with the appropriate hyperangular, angular, spin 
and isospin wave functions f34|| . In particular the nucleon wave function is written as 

\N, J p = 1/2+) =^oo(0^°5] ) (^)X5(l/2;l/2) , (50) 

where xs(l/2; 1/2) is the symmetric SU(6) combination of spin and isospin wave functions of the three quarks, and we 
have explicitly indicated the symmetry required for the hypershperical wave function and the total angular momentum 
L = l p + l\: symmetry] • -^ or an y gi ven SL/(6)-configuration only ipv,f(0 will depend on the potential and we can 
write general results for all the sum r ules when H = T + V^(0 §|. 



1. sum rules for the electric dipole excitations 



i) m (D' z ): 

m (D' z ) = (J2 eUf) + {J2 e i e j z' i z^) = 

i i^j 

= \e 2 J dU 7 \^oo(0\ 2 = le 2 (rl) ch (51) 

both for protons and neutrons. The non-energy- weighted sum rule is therefore proportional to the mean square mass 
radius of the nucleon, which, at least as long as the nucleon wave function belongs to a SU(6) multiplet, is also the 
proton charge radius. 

ii) mi(-D') = e 2 /3m remains unmodified for all the velocity independent potentials as already discussed (cfr. Eq. (|32])) 
and does not distinguish neutron and proton. Its model independence is based on the property [V( 3 )(£), D' z ] — 

iii) m 2 (D' z ) 

m 2 {D' z ) = < i( ei + e 2 2e 3 ) 2 V\ V\ + ( ei - e 2 f V* V* + 

+ -j= (ei + e 2 - 2e 3 ) (d - e 2 ) V z p V z x ) = 

= ^ < T > = + f *»<0 • («> 

77i2 does not depend on the interquark interaction explicitly, however is proportional to the mean kinetic energy of 
the system and measures the presence of high momentum components in the nucleon wave functions. If the potential 
is highly confining the high momentum components are larger and the sum increases. 



iv) m 3 (D' z ) 



m 3 (D' z ) = ^{\ (ex + e 2 - 2e 3 ) 2 (V* A V{ (£)) + 
2m z o 

+ ^= (ei + e 2 - 2e 3 ) (ei - e 2 ) (VJ V* V&> (£))) = 



18m 2 J " \ d£ 2 C d Z 

The three-times energy-weighted sum is particularly interesting because it depends crucially on the potential model 
and on its derivatives. 







2. sum rules for magnetic dipole excitations 



Magnetic excitations induced by the dipole operator (j^) are easily understood in the case of spin-independent 
interactions. In fact the magnetic operator commutes with the nucleon Hamiltonian [T + V' 3 )£), fiz] — and the total 
strenght is concetrated at a value of the excitation energy equal to the ground state and one gets: 

i) moinz) 



m (fi z ) = 



(2 m) 2 
\2m) 



e, e at a]} + e 2 ) - | e, a 



z\ |2 
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1 - 1 



Mo 



(54) 



The result (|5J) is valid for both protons and neutrons (in the neutron case the indidual contributions in parenthesis 
become 2/3 + 2/3 — 4/9), and we have defined /iq = e/2m. 

ii) higher sum rules: 



mi(n z ) = m 2 (p z ) = m 3 (/i 2 ) = 



(55) 



vanish because we did not consider the effects due to hyperfine interaction. In the following we discuss a variety of 
hypercentral potential investigated in ref. |^7j , and to this hand we discuss first the role of the spin-spin interaction 
terms. 



D. The role of the hyperfine interaction 



The authors of the ref. |27| ] include a perturbative Hyperfine interaction in the three-body Hamiltonian in order 
to reproduce the correct iV-A mass splitting. Such contribution is seen as the non-relativistic reduction of the one- 
gluon-exchange interaction j37| . We restrict the Hyperfine contribution to the dominant spin-spin zero-range term 



T/ (2) / s 2 as 8 7r 
^Hy P ( r i2) = — Si • S 2 <5(r 12 ) 



3 m 2 3 

neglecting the (small) tensor contribution. The effective coupling constant as is fixed by the N-A mass splitting 

V2 a s 8 



(56) 



M A -M N = ^^* fdU 2 \MO\' 
6 m z tt J 



(57) 



The additional contributions due to the Hyperfine potential (56) prevents the magnetic dipole sums to vanish, and 
modifies the energy-weighted electric dipole sums. 



1. sum rules for the elctric 



e excitations 



(2) 

zi "Hyp 



i) mi(D' z ) and m,2(D' z ) remain unmodified since D',, V T 

ii) m^{D' z ) gains an additional contribution: 

\f2ir as 



™ S (D' Z )\ 



Hyp 



((e 1 -e 2 ) 2 S 1 -S 2 [V 2 p 6(p)}) 
- - (Ma - M N ) J W¥¥ ^ m 



9 TO 4 

1 1 2 
3 m 



(58) 



which has to be added incoherently to the Eq. (p3|). In Eq. (^8|) i^ooiO indicates the total derivative of the wave 
function and the relation (|57|) has been used J36|. 
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2. sum rules for magnetic dipole excitations 



The crucial role of the Hyperhne interaction for the magnetic energy- weighted sum rules has been already pointed 
out in the previous sections by stressing the fact that these sum rules simply vanish if the spin-spin term J56| ) is 
neglected. Hyperfine potential is, however, a two-body operator and the commutators with the magnetic dipole 
moment (||) are not as simple as the electric case. For simplicity we will restrict to the linear energy-weighted sum 
rule which, differently from mi(D' z ), is not model independent and embodies the interesting aspects of the interaction- 
dependence we recognize in my,{D' z ). 

As long as the eigenstates of the hypercentral potential are assumed to belong to a specific SU(6) multiplet, the 
magnetic operator can excite one state only: the A33 resonance. Within such simple N-A excitation model the 
magnetic susceptibility is given by 

_ ](A|^iiV)i 2 _ 16 2 1 4 a 

M A -M N - T^Ma-A/*- 8 - 7 ' 10 fm • (59) 

The sum rule approach cannot add any information to this simple scheme as it can be demonstrated calculating the 
first moments of the magnetic strength distribution: 

i) m (fi z ) has the form (|54|). 

ii) mi(p z ) 

mM = ~ ((ei - e 2 f (S 1 ■ S 2 - S{ Sf) ^ ^ 6(r 12 )) 

= |**o ( m a - M N ) , (60) 

where Eq. ( |57j) has been used. 

By using the bounds involving mi , too and E±q = M A — Mm one obtains 

,2/ 



o "W z ) _ 16 2 1 s a s o m o(^) _ 16 2 1 r l 

~Y^M A -M N - (jM - 2 ^r~Y fl0 M A -M N ' 



and the two limits coincides with the result (^9|) as expected. 

An estimate of (3m different from the previous result comes from SU (6) breaking components in the nucleon wave 
function. In the next section we discuss a potential model whose nucleon wave function results in a superposition of 
different SU(6) configurations. 



E. Two-body + three-body interactions 

The pathologic aspect of non-relativistic quark models involving confining potentials plus (at least part of ) the 
one-gluon-exchange interaction remains the impossibility to predict the masses of the first (negative parity) excitations 
in the baryonic octet and decuplet correctly. The authors of ref. ||| proposed to add a phenomenological three-quark 
force getting a correct spectrum up to excitation energy of 0.7 GeV. The analysis of the mesonic qq and baryonic qqq 
spe ctrum leads to a two-body potential containing confining, Coulomb-like and spin-spin qq interaction of the form 

,(2) _ 1 / ^ + !^2 +4 ^£_ exp(-ri 2 /r ) 

2 \ ri2 a 2 to 2 rl r 12 



V i-i = 9 + + 4 ^2 ^7 Si ■ S 2 - £> ) , (62) 



where the Yukawa form of the spin-spin term replaces the delta contact interaction of the OGE potential ( |56|) to avoid 
an unbounded spectrum when the Scrodinger equation is solved. 
The three-body term suggested by Cano et al. in ref. J3S" 

vfU = V o ex P I - U TJ I = y o ex P (-3£7Ao) ( M • 




11 



(3) 

is identified as Vjjj by the same authors and assumes a simple form within the hyperspherical formalism. The 
parameters of the interaction can be found in table 1 of their paper and in the following we discuss the results of such 
quite recent version of the QCD inspired potential including two- and thre-body contributions. 

The solutions of the Schrodinger equation can be expanded on the hyperspherical basis and this is the actual 
procedure followed in ref. |B2| . In practical calculations only two terms have been retained and the nucleon wave 
function can be written: 

\N, J p = 1/2+) = Y^ff X[s] (1/2; 1/2) + 

+ tf 3 (0 -L [y^ X[MS] (l/2; 1/2) + Y$$ A] X[MA](l/2; 1/2)] , (64) 

where the second term of the expansion is clearly an S'f7(6)-breaking contribution, and X[MS](l/2; 1/2), 
(X[ma](1/2; 1/2)) are the mixed-symmetric (antisymmetric) SU(6) combinations of spin and isospin wave functions 
of the three quarks. Owing to the 5t/(6)-breaking term in the nucleon wave function, some results of the pure hy- 
perradial potentials are modified. For istance the charge root mean square radius of the neutron is not vanishing and 
the values depend on the contributions coming from the three-body part of the potential. 



1. sum rules for electric dipole excitations 



i) m (D' z ): 



m (D' z ) = le 2 y d^ 7 [|*i(£)| 2 + |* 3 (£)| 2 ] -^=J d^^MOMO 



\e 2 (r 2 ) ch + ^e 2 (r*) ch , 



(65) 



both for protons and neutrons. The non-energy-weighted sum rule, therefore, is no longer proportional to the charge 
mean square radius of the proton when SU(6) breaking component in the nucleon wave function are taken into 
account. 

ii) mi(D' z ) — e 2 /3m again remains unmodified because [V^_ (£) + (£ ), D' z ] — as already discussed and does 
not distinguish neutron and proton. 

iii) m2{D' z ) is no longer proportional to the mean kinetic energy and contains an additional term: 



, ,s 2e 2 , v e 2 1 
=9^ + 9^71 



d 2 1 d 



(66) 



iv) m 3 (D' z ) 

Both V q _ q and V q _ q _ q contribute to the sum rule which can be written m 3 {p z ) = m i{D' z )\ q _ q 



where 



"*ra g - 9 = ^ <(ei-e a ) 



2 V 2 



«C . f 12 . . k ct exp(-r 12 /r Q ) 

1 7T + 4 — o 9 S l ' S 2 

r 12 a z Tfi r 12 



m 3 (D' z )\ q _ q _ q , 



(67) 



and the three-body m 3 {D' z )\ is again given by Eq. (||) with the obvious replacement F (3) (0) ~> Vg- q -q(0) 



2. sum rules for magnetic excitations 

The role of SU (6) breaking components in the nucleon wave function ( |6~i| ) is particularly relevant for the magnetic 
excitations because the operator \x z can now mix different SU(6) configurations. We discuss the simple bounds which 
involve m , mi and only. One gets 

i) m (n z ) 



mo(ji z ) = Mo 



9 Pl 



ft + 3ft 



(68) 
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where P a — f d£ £ 5 \ip a (0 1 2 an d one recovers the result ( |54"| ) in the limit Pi = 1 and P3 = 0. In the neutron case the 
sum rule reads 



1 ( Px + ift) - §p? 



(69) 



However, despite the different form of the sum rules the numerical values are quite close each other and do not produce 
differences in the estimates of the /3m for protons and neutrons. 

ii) mi(p z ) 



mi(p z ) 



12 



(( ei - e 2 ) 2 (Si • S 2 - SI S z 2 



1 K a exp(-ri 2 /r ), 



777 z 4 777 



(70) 



IV. NUMERICAL RESULTS AND DISCUSSION 

We discuss, first, numerical results for a variety of hyperradial potentials introduced by Ferraris et al. f27j] , namely 

(71) 
(72) 

(73) 
(74) 

(75) 

whose parameters are summarized, for convenience, in table [l| The potentials V\ - V3 have been fitted on the baryon 
mass spectrum, the hypercoulomb parametrization (|74|) has been considered as check of numerical calculations, while 



v} 3 \o = 


T 

~Y 


vhi 




v 2 {3) (0 = 


T 

~Y 




b 

~e 


v 3 (3) (0 = 


T 

~Y 




b 


vf\0 = 


T 








T 

~Y 


b 

h F 





clogC 



the version ( |75| ) has been introduced in ref. [ 
square radius of the proton. In the same tab 



27] because it reproduces the electric dipole form factor and root mean 
e the predictions of the mass radius are also presented. 



TABLE I. Parameters of the potentials ([n]) - (JtH). In the last coulomn the predicted mass radii. 



potentials 


T 


ki 


b 




<r 2 ) 




M 


[fin" 2 ] 


[fm] 


[fa- 1 ] 


[fm 2 ] 




4.59 


1.61 






(0.516) 2 




2.50 


1.14 


-0.80 




(0.483) 2 


v 3 (3) 


2.50 


1.14 


-0.80 


0.1 


(0.475) 2 


n (3) 


6.39 








(0.462) 2 




1.78 




-0.78 




(0.88) 2 
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A. Electric dipole polarizability 



1. results of three-body force models (TBM) 

The electric dipole sum rules are shown in table || for the hyperradial interactions ( |7l"| ) - (fT5|) , and the corresponding 
lower and upper limits in tables II], and IV. The smaller is the predicted radius of the system, the smaller is the 
non-energy- weighted sum and the larger the mi moment. In fact mo(D' z ) is proportional to the proton charge radius 
and m2{D' z ) to the mean kinetic energy which is larger in smaller systems because of the indetermination principle. 
The role of the hyperfine interaction is clearly seen in m^{D' z ) (cfr. table ||). A large part of the sum comes from 



the spin-spin interaction and the lower limit 2mf/m3 can become significantly small (cfr. table III) loosing a direct 
connection with reasonable values of the polarizability. One can conclude that the charge deformation induced by the 
external electric field on the quark distribution is far from beeing approximated by a rigid oscillation of u and d charge 
densities (opposite in phase) as assumed in the relations ([f2| ) , and ([h|) . The inclusion of hyperfme contributions to 
the m.3 sum lowers the polarizability by few percent enlarging the difference between the upper and lower limits. 

TABLE II. Electric dipole sum rule values for the potential models (|7l|)-(|75|) (cfr. also table |). For the 7773 sum the 
contributions coming from the hyperradial part of the potential is shown, the total result contains, in addition, the Hyperfine 
contribution coming from the interaction term (|56|). 



potentials 


m (D' z ) 
[mi 2 ] 


mi(D' z ) 
[fm] 


m 2 (D' z ) 

M 


m s(£>z)lv(3) 
[fm" 1 ] 


7713 (-Dz) tot 

[fm- 1 ] 




0.0006 


0.0015 


0.0040 


0.0137 


0.0278 




0.0006 


0.0015 


0.0056 


0.0760 


0.2085 




0.0005 


0.0015 


0.0058 


0.0799 


0.2149 




0.0005 


0.0015 


0.0053 


0.0274 


0.0498 


V(3) 


0.0019 


0.0015 


0.0018 


0.0098 


0.0602 



TABLE III. Lower bounds to the electric polarizability as predicted by the potential models (|7l|)-(|75]). In parenthesis the 
numerical results obtained neglecting the Hyperfine contribution in 777,3 (cfr. table 0). The polarizability is expressed in 10 
fm 3 . 



potentials 


>"?(£>i) 


m f,(K) 

* mi(Ot) 


m l( D ',) 
z mi (i5C)(i-A/r) 




1.69 (3.44) 


5.48 


5.50 (5.60) 


vf> 


0.23 (0.62) 


4.19 


4.22 (4.29) 


vf> 


0.22 (0.59) 


3.94 


3.97 (4.04) 


vf> 


0.94 (1.72) 


3.51 


3.55 (3.65) 


vf> 


0.78 (4.79) 


45.7 


45.9 (47.1) 
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Looking at the comparison between the simple lower bound 



(76) 



shown in table III and the results of table fv| one can see that the approximation ((76|) is rather good for all the potential 
models suggesting that Eq. (Q) with F = D' z represents a more reliable nucleon charge deformation induced by the 
electric field. 

TABLE IV. Upper bounds to the electric polarizability as predicted by the potential models (0)-(^). The polarizability is 
expressed in 10 -4 fm 3 . 



potentials 


m\{D z) m.Q(D' z) \— • 
rn 2 (D f z) 


1 m (D') 


2 2#iA 


2y/mi(D' z )m- 3 (D' z ) 




5.65 


5.86 


5.71 


5.88 




4.53 


5.08 


4.78 


5.04 




4.27 


4.79 


4.51 


4.75 


n (3) 


3.64 


4.09 


3.82 


3.95 




46.5 


66.2 


58.1 


56.0 



TABLE V. Bounds to the electic polarizability as predicted by the potential models ([7l])-(|7^). The results obtained neglecting 
Hyperfine contributions are also shown. The polarizability is expressed in 10~ 4 fm 3 . 



potentials 


V (3) (0 [only] 


v< 8 >(0 + ^ 




a!; 11 = 5.65 ± 0.05 


aP' n = 5.60 ±0.10 




a P' n = 4.54 ± 0.25 


a^ n = 4.50 ± 0.28 


v 3 (3) 


aP' n = 4.27 ± 0.23 


aP' n = 4.24 ± 0.27 


vf> 


a% n = 3.74 ±0.09 


a% n = 3.69 ±0.14 


vf> 


«P' n = 52.6 ± 5.50 


ftP' 11 = 52.9 ±6.10 



2. results of Two-body + three-body force model 



In the following we present results of two sets of potentials (cfr. section |IIE|): i) Vj: includes two-body interaction 
only (see Eq.(^32|)); Vjji'- includes a three-body term also as in Eq. (|6^). The numerics is summarized in tables [vj ■ 



[X, where also the effects of the hyperfine interaction and three-body terms are emphasized. In particular the role of 



three-body interaction i s qu ite important. Neglecting that contribution to my,{D' z ) yields to a negative value of the 



lower bound (cfr. table VII), while the hyperfine interaction plays a role similar to that described in relation with 
the hyperradial models. The overall impression is that, differently from the hyperradial models, the two body models 
have a quite small radius and therefore a smaller electric polarizability. Adding the three-body term does not improve 
the situation and the resulting polarizability is even smaller. 

TABLE VI. Electric dipole sum rule values for the potential models Vjr, Vni (|62|) - (j63|) . For the 7713 sum the contributions 
arising from the various potential terms of Eq. (^) are shown separately; the last value refers to the three-body term and it 
vanishes for the Vi version. 



potentials 


m (D' z ) 


mi(D' z ) 


m 2 (D; 




m 3 (^) tot 




[fm 2 ] 


[fm] 


M 




[fm- 1 ] 


Vi 


0.0005 


0.0014 


0.0044 


0.0168 


= 0.0022 ± 0.0073 + 0.0072 ± zero 


Vni 


0.0003 


0.0015 


0.0197 


0.7441 


= 0.1214 ± 0.0040 + 0.0524 ± 0.5662 
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TABLE VII. Lower bounds to the electric polarizability as predicted by the potential models Vi, Vm (|6j)-(|63j). The results 
in brakets () for the version I are obtained neglecting the Hyperfine contributions; those ones in [] for the version III are 
obtained excluding the three-body contributions. The polarizability is expressed in 1CP 4 fm 3 . 

potentials 2 2 "Aff-* 

_ m 3 (D' r ) mi(g;) ">i(D',)(l-A/r) 

Vi 2.42 (4.25) 3.27 3.30 (3.24) 

Vm 0.06 [0.25] 0.97 1.15 [-0.061] 



Another interesting conclusion can be drawn for the approximated expression (|4q). Its values are still not far from 
the exact results, but the inclusion of three-boy forces worses the situation. One has to emphasize, however, that 
the bound 2m§/mi cannot be expressed in terms of the mean charge radius of the proton because of the SU(6) 
breaking components present in the wave function ( |64[ ) and the second identity ( |76| ) is no longer valid. If one assumes 
Qig n = |e 2 (3m) (rp) c h one should have = 4.37 10 -4 fm 3 for Vi and = 1.12 10~ 4 fm 3 for Vm, values rather 
far from the exact result of table [X: the ST/ (6) breaking components produce relevant effects on the polarizability 
which are not included in the most simple estimates. In particular the bounds which involves higher sum rules (7713 
and m%) are quite sensitive to the tuning effects due to the ^3 components and if one puts ^3 = in the calculations 
one cannot satisfy the correct inequalities obtaining, for the Vi potential model, a lower bound (cxe = 3.82 10~ 3 fm 3 ) 
larger than the upper bound (ajg = 3.73 10~ 3 fm 3 ). 
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TABLE VIII. Upper bounds to the electric polarizability as predicted by the potential models Vi, Vm (p2[)-(|63[). The 
polarizability is expressed in 10~ 4 fm 3 . 

potentials 2 ZSgil 2 A 

£m Eia 

Vi 3.58 3.38 

Vm 1.87 1.64 



TABLE IX. Upper and lower bou nds to the nucleon polarizability as predicted by the potential models Vi, Vm (|6^)-(|63|) 

The polarizability is expressed in 10~ 4 



Notations for the brakets as in table VII 



fm 13 



potentials 



Vi 
Vm 



a'g' 11 = 3.34 ± 0.04 
a p ;° = 1.39 ± 0.24 



3.31 ± 0.07) 



[ag n = 1.42 ± 0.45] 



B. Magnetic susceptibility 

Paramagnetic susceptibility of hyperradial models does not differ from the results of Eqs. (|5^) and ( |6l"| ) because the 
Hamiltonian (^) commutes with fi x ■ On the contrary the SU (6) breaking model with three-body forces has a much 
more complex spin structure which is evident in Eqs. (|68[), and (|70|). The results are shown in tables [x| - XIII. 



TABLE X. Values of the magnetic sum rules for the potential models Vi, Vm (|62[)-(| 
potentials ma(^z) 

[M H 

Vi 5.60 10~ 4 0.0011 

Vm 6.17 10~ 4 0.0015 



TABLE XL Lower bounds to the paramagnetic susceptibility as predicted by the potential models Vi, Vm (p2[)-(|63|). The 
susceptibility is expressed in 10 -4 fm 3 . 

potentials 2 m "^'j 

Vi 5.89 
V/ii 5.23 



Since we are calculating the simplest bounds only, the uncertainties are larger than those ones for the polarizability 
(13% and 23% for Vj and Vm respectively), however the reduction with respect the N-A approximation is quite 
relevant (roughly —30%) and is due to inclusion of the 7 = 2 component in the hyperspherical expansion (|&i|). Again 
the results are identical for protons and neutrons. 

TABLE XII. Upper bounds to the paramagnetic susceptibility as predicted by the potential models Vi, Vm (|6^)-(^). The 
susceptibility is expressed in 10 -4 fm 3 . 

potentials 2 m ° ( fl M / ) 

fllD 

Vi 7.55 
Vm 8.31 
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TABLE XIII. Upper and lower bounds to the paramagnetic susceptibility as predicted by the potential models Vj, Vm 
32|)-(B3). The susceptibility is expressed in 10~ 4 fm 3 . 



potentials 


Vi 


oP,n _ 
t J M " 


= 6.72 ±0.83 






= 6.77 ± 1.55 



C. Comparison with experiments and conclusions 



In tables |XIV . XV all the results are summarized and the retardation corrections included for both protons and 
neutrons. The calculation of the additional terms, within the quark models we are discussing, is rather straightforward 
because they are basically related to the electric dipole sums already calculated in the previous sections and to the 
charge proton and neutron radii. In particular (cfr. Eqs. (|J)) and (^)) 



e 

3M 



p,n 



'ch 



(77) 



e 

2M 



/ch 



(78) 



A ldia 6m 



(79) 



e 

6m 



)ch + -j (rj 



ch 



(80) 



where also the SU(6) breaking effects have been included. 
Let us discuss first the electric polarizability. 

The values closest to the experimental data are obtained for the hyperradial potential models, where the three- 

(3) 

body contributions are included in the structure of the interaction. The peculiar potential V5 introduced because 
it reproduces the electric form factor of the proton, shows the classical limitation of the quark models: a huge 
polarizability when one ask to the quark degrees of freedom to cover all the charge spatial distribution replacing 
also the effects of the meson cloud. This is a well known problem and it can be easily understood looking at the 
approximate relation (|35|): by replacing the charge root mean square radius with its experimental value one gets 
cle ~ 42 10~ 4 fm 3 . On the contrary if one wants to reproduce the excitation energy of the dipole states a>h.o. ~ 600 
MeV, and one obtains ~ 3.3 10~ 4 fm 3 ! In order to overcome this contraddiction we want to investigate better the 
predictions of the quark models which are able to reproduce the spectrum of baryons and look at the quantitative 
amount of meson contribution one needs to fill the gap between quark model contribution and the experimental data. 



TABLE XIV. Retardation corrections to the electric polarizability 
the hyperradial potential models (|7^) - (fr3) and the models (|6^)-(|63| 



and magnetic susceptibility of the proton calculated for 
. The units are 10~ 4 fm 3 . 



potentials 






A/3 P 




Pit 


V< 3) 


1.36 


6.96 


-2.04 


-2.04 


4.61 


vf> 


1.19 


5.70 


-1.79 


-1.79 


5.13 


vt 


1.16 


5.40 


-1.73 


-1.73 


5.23 


v 4 (3) 


1.09 


4.78 


-1.64 


-1.64 


5.24 


vf> 


3.93 


56.8 


-5.90 


-5.90 


-3.10 


Vi 


1.22 


4.56 


-1.41 


-1.60 


3.71 


Vm 


0.61 


2.00 


-0.83 


-0.87 


5.06 
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In the case of hyperradial TBM potentials let us take the results coming from the linear plus confining " Coulomb" 
(i.e. the interaction ) as the typical example. Which contribution should one need from mesons to reproduce the 
experimetal polarizabilities of neutron and proton? Including the meson contributions into the retardation corrections 
( |77| ) is rather simple by taking the experimental value of the proton charge radius: one obtains the prediction for the 
dynamic polarizability aP E = (5.60 + 3.78) 10~ 4 = 9.38 10~ 4 fm 3 which is claming for a contribution of the mesons 
in the range of ss 2.7 10 -4 fm 3 . For the neutron the static value of the polarizability is known experimentally (cfr. 
section ffl) and the meson contribution is found to be much larger ps 7.3 10 fm 3 in agreement with the calculation of 
ref. [|15|| where the difference between the neutron and proton static polarizability is found 5 10 -4 fm 3 and ascribed 
to mesonic degrees of freedom. The net result is that the contribution of the quarks to the proton polarizability is 
not small and comparable with the one coming from the meson cloud. This conclusion is basically valid for all the 
TBM Vi — fitted on the excitation spectrum and it favours the QCD "inspired" form — t/£ + fc; £, i.e. the v} 3 ^ 
potential model. On the contrary the potential models where the three-body forces are included phenomenologically 
to reproduce the position of the Roper resonance, predict quite low values of the static polarizability therefore asking 
for a much larger meson contribution j42] | . 

The results on the magnetic susceptibilities 

The predicted paramagnetic susceptibility is the same for all the TBM potentials /3?' f n I Dara = 8.7 10~ 4 fm 3 . In 
particular for the proton the inclusion of meson spatial distribution on the correction ( |78|) and on the diamagnetic 
contribution ([79]) yields to J3 P M = (8.7 — 9.0) 10 -4 fm 3 and therefore to a mesonic contribution 2.4 ± 1.3 10~ 4 fm 3 , a 
result consistent with the analysis of ref. where the Ml pion photoproduction has been investigated. 



TABLE XV. Retardation corrections to the electric polarizability and magnetic susceptibility of the neutron calculated for 
the hyperradial potential models (173) - ( fr5| ) and the models (S)-(B3|). The units are 10 -4 fm 3 . 



potentials 


Aa n E 




Af3 n 


fldta 


0u 







5.60 


-2.04 


-1.36 


5.29 







4.50 


-1.78 


-1.19 


5.72 







4.24 


-1.73 


-1.16 


5.81 


n (3) 





3.69 


.1.64 


-1.09 


5.97 







52.9 


-5.90 


-3.93 


-1.16 


Vi 


-0.10 


3.24 


-1.41 


-0.99 


4.32 


Vm 


-0.02 


1.37 


-0.83 


-0.57 


5.37 
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The neutron susceptibility has been experimentally derived by the knowledge of the static electric polarizability, 
adding retardation corrections and using the dispersion relation (||), the precision is rather limited. The mesonic 
contribution can be extracted from /3" f = (8.7 — 6.8) fm 3 and it results to be 1.9 ± 3.6 10~ 4 fm 3 . 

The potential models Vj - Via give analogous results. The remarkable difference is that their complicate spin 
structure can be easily seen looking at the mean excitation energy ml(/i z )/mo(/i z ) ~ 370 MeV and 467 MeV respec- 
tively which is larger than the Ma — Mm indicating that the SU (6)-breaking mechanis m int roduces a larger excitation 



spectrum: an effect which results also in the width of the predictions for /3m in table XIII 
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